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High resolution 21 cm observations of the Ursa Major cirrus revealed highly filamentary struc- 
tures down to the 0.03 pc resolution. These filaments, still present in the line centroid map, show 
( multi-Gaussian components and seem to be associated with high vorticity regions. Probability 

density functions of line centroid increments and structure functions were computed on the line 
centroid field, providing strong evidences for the presence of turbulence in the atomic gas. 
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Many statistical studies of the density and velocity structure of dense interstellar mat- 
, ter have been done on molecular clouds where turbulence is seen as a significant support 

against gravitational collapse that leads to star formation. Less attention has been de- 
voted to turbulence in the Galactic atomic gas (HI). The cold atomic component (T ~ 100 
K, n ~ 100 cm -3 ), alike molecular gas, is characterized by multiscale self-similar struc- 
tures and non-thermal linewidths. 
' A detailed and quantitative study of the turbulence and kinematics of HI clouds has 

never been done. Here we present a preliminary analysis of this kind based on high 
resolution 21 cm observations of an HI cloud located in the Ursa Major constellation. 
To characterize the turbulent state of the atomic gas, a statistical analysis of the line 
centroid field has been done. We have computed probability density functions of line 
. . \ centroid increments and structure functions. 
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^ ' 2. HI Observations 

H 

C3 ; The Ursa Major cirrus (a(2000) = 9ft36m, (5(2000) = 70°20') has been observed with 

the Penticton interferometer. Half of these data has already been published by Joncas et 
al. (1992). A total of 110500 spectra were taken on a 3.75° x 2.58° field (1' resolution) 
with a spectral resolution of Av = 0.412 km s _1 . This cirrus is part of a large loop of 
gas, possibly blown by stellar winds and/or supernovae explosions, known as the North 



Celestial Loop. A distance of 100 pc is adopted ( de Vries et al. 1987 ). 

The integrated intensity map, proportional to the column density, is shown in figure 1; 
many filamentary structures are apparent. The filaments stand out in individual channel 
maps and most of them can be followed over several channels. The strongest filament (a 
on figure 1) crosses most of the field in an east-west direction. There are also two fainter 
filamentary structures (b and c) that seem to merge (in spatial and velocity space) in 
the brighter part of filament a. A closer look at these filaments in individual channel 
maps reveals smaller elongated substructures, parallel to the main axis. 

We performed a multi-Gaussian fit on each spectrum and found that most spectra are 
fitted by one or a combination of two Gaussians of typical velocity dispersion a = 2.8 
km s _1 . The smallest velocity dispersion measured (a ~ 1.3 km s -1 ) is located in the 
northern part of filament c. If this width were purely thermal it would correspond to a 
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Figure 2. Map of line centroids for spectra with a signal-to-noise ratio greater than 4. 
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Figure 3. Normalized Probability Density Function of velocity increments (lag of 3 and 50 
pixels) computed from the line centroid field of figure 2. 

kinetic temperature of 200 K. The spectra along filament a are composed of two Gaussian 
components (see central spectrum on figure 1). These two components are separated by 
Ai> = 4.0 km s _1 all along the filament except in the eastern end where they merge. At 
this position we observe a clear velocity gradient (10.0 km s _1 per pc) across filament a. 
These two observational facts may be interpreted as a rotation in the plane of the sky, 
along the major axis of the filament. Assuming the filament is a cylinder, its radius is 
R = 0.2 pc. Then, the turnover time (rj = 2R/Av), calculated from both behaviours 
(component separation and velocity gradient) is r t = 1 x 10 5 year. 



3. Statistical analysis of the line centroid map 

Figure 2 displays the map of the line centroid (vo — J v Tb(v) dvj f Tb(v) dv). On large 
scales, the line centroid field is characterized by a global east-west velocity gradient. On 
small scales, we note a significant difference between the eastern part (particularly near 
a = 9h47m, S = 71° 15') which shows velocity fluctuations on scales of a few arcminutes 
and the western part which is very smooth at those scales. The brighter part of filament 
a as well as filament b and c clearly stand out as coherent structures in the centroid map. 

3.1. Probability Density Function of velocity increments 

Probability density functions (PDF) of velocity increments and derivatives found in ex- 
periments and numerical simulations of turbulent flows usually show clear non-Gaussian 
wings. This non-Gaussian behaviour is attributed in part to coherent vortex structures 
( [Breenivasan fc Antonia 1997 ). In astrophysical context, non-Gaussian wings were ob- 



served by Falgarone & Phillips (1990) in molecular line profiles of selected molecular 
clouds. Recently, Lis et al. (1996) showed that PDFs of line centroid velocities incre- 
ments (i.e. PDF of Svq = vo(r) — vq(t + Sr)) computed from numerical simulations of 
compressible turbulence at high Reynolds number have non-Gaussian wings. 

Using the line centroid map shown in figure 2 we computed PDF of velocity increments 
for 5 and 50 pixels lags. As shown in figure 3, non-Gaussian wings are clearly present in 
the Ursa Major cirrus. Furthermore, exponential wings are stronger at small lag. This 
is a behaviour typical of intermittency also observed in three dimensional simulations 



of turbulence and in laboratory experiments ( |5reenivasan fc Antonia 1997 ). The points 
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Figure 4. Second and third order structure function of the line centroid field. The structure 
functions were computed on the west and east side of the field 



that populate the non-Gaussian wings are essentially locate d along filaments a, b and 



c, and in the north-eastern part of the field. In their analysis Lis et al. (1996) showed a 



clear spatial correlation between the non-Gaussian events and regions of high vorticity. 

3.2. Structure Functions 

Structure functions (SF) of various orders (< Au n >— -j^Yl ^ v q) are extensively used to 
deduce physical parameters of turbulent flows (Ansclmet et al. 1984). Second and third 
order SF were computed from the line centroid map. The large number of velocity points 
used here allows us to separate spatially our sample in two equal parts (east and west). 
We found a significant difference between the eastern and western SFs. 

In the inertial sub-range of a Kolmogorov type turbulence, the exponent of the second 
and third order SF should be respectively 2/3 and 1. Here, as seen in figure 4, the second 
order SFs are very smooth and show various exponents (from n2 = 0.25 to ri2 = 0.89) 
on different scales and locations. The second order SFs drop near 1.5 pc; this may be 
caused by the poor statistics we have at large scale. The third order SFs are also very 
smooth. The exponents range from n% — 0.96 to 713 = 2.20. 

From 0.4 pc to 1.5 pc, the eastern part of the line centroid field is characterized by 
slopes of ri2 = 0.77 and rc.3 = 0.96, and seems to meet the Kolmogorov type description 
for which the third order SF has an exact solution: 



< Ai 



>= --r <e> 
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where < e > is the mean energy transfer rate in the inertial range. This rate can be 
computed for the eastern part of the field: < e >= 1 x 10 -2 erg s _1 g _1 . On a small 
lag range (from 0.4 pc to 0.8 pc) we also observed an exponent of 1 in the western third 
order SF, giving < e >= 0.8 x 10~ 4 which is significantly smaller than the eastern value. 
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4. Discussion 

The previous statistical analysis supports the presence of turbulence in the Ursa Major 
cirrus. The various observational results shown here (the presence of clearly non-Gaussian 
wings in PDFs of line centroid velocity increments and the slopes of the structure func- 
tions) are recognized as signatures of turbulence in terrestrial experiments or numerical 
simulations. On the other hand, astrophysical observations lack three dimensions in 
phase space; any detailed comparison with numerical simulations or laboratory results 
must be done with great caution. Here the 3D geometry of the cloud is unknown, it is 
thus premature to interpret further the various SFs' exponents. 

Numerical simulations and laboratory experiments of turbulent flows are also charac- 
terized by filamentary vortical structures and dissipative sheets around them 
( [Breenivasan fc Antonia 1997 ). These highly localized dissipative regions may play an 
important role in the thermodynamics and chemical evolution of the interstellar medium 
( Falgarone fc Puget 1995 ). Identifying such regions would be of considerable astrophysi- 
cal interest. The present study of high resolution 21 cm observations reveals filamentary 
HI structures characterized by doubled-peak spectra and minor axis velocity gradient. 
These features may be interpreted as a rotation around the major axis of the filament. 
Furthermore, the main filaments observed here are associated with high vorticity zones 
traced by the PDF of velocity increments. 

The filaments we have observed in the Ursa cloud cannot be the coherent vortices which 
develop at scales close to the dissipation scale of turbulence because the linear resolution 
of our observations (0.03 pc) is more than 100 times larger than the dissipation scale in 
diffuse clouds (see Falgarone, these proceedings). Instead, we may be observing a larger 
structure, made of very small vortices, braided together, as the result of their interaction 
and merging. This merging of small vortices into larger vortices has been followed in the 
numerical simulations of Vincent and Meneguzzi (1994) and Porter et al. (1994). This 
would explain why the velocity pattern is more complex than a simple rotation. 
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